Spin injection and electric field effect in degenerate semiconductors 
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We analyze spin-transport in semiconductors in the regime characterized by T ~ Tf (intermedi- 
ate to degenerate), where Tf is the Fermi temperature. Such a regime is of great importance since 
it includes the lightly doped semiconductor structures used in most experiments; we demonstrate 
that, at the same time, it corresponds to the regime in which carrier-carrier interactions assume a 
relevant role. Starting from a general formulation of the drift-diffusion equations, which includes 
many-body correlation effects, we perform detailed calculations of the spin injection characteristics 
of various heterostructures, and analyze the combined effects of carrier density variation, applied 
electric field and Coulomb interaction. We show the existence of a degenerate regime, peculiar to 
semiconductors, which strongly differs, as spin-transport is concerned, from the degenerate regime 
of metals. 

PACS numbers: 72.25.-b, 72.10.-d, 72.25.Dc 
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I. INTRODUCTION 

Despite many efforts and some substantial progress, 
major difficulties persist in the practical implementation 
of semiconductor spintronics devices. One of the funda- 
mental problems inherent the electrical injection of spin 
currents into semiconductors is the so-called conductiv- 
ity mismatch^: the large conductivity difference between 
metals and semiconductors induces a spin accumulation 
at the interface which in turn drastically reduces spin 
injection into the semiconductor. 

Recently Yu and Flatte (YF)— have pointed out the 
possibility of improving the efficiency of spin injection 
through the application of a strong, but experimentally 
accessible, electric field. At variance with metals, in 
semiconductors, this field remains partially unscreened so 
that it is necessary to explicitly introduce a drift-related 
term in the spin-transport equations. 

The YF analysis of the semiconductor high-field regime 
neglects electron-electron interaction effects - aside from 
those responsible for macroscopic charge neutrality - 
and, perhaps more importantly, focuses on the high-field 
non-degenerate regime (see Eq. l(T]l). especially insofar as 
the analysis of heterostructures are concerned. 

In the present paper we extend and refine their analysis 
in several directions. Starting from a more general for- 
mulation of the drift-diffusion equations, which includes 
many-body correlation effects, we perform detailed cal- 
culations of the spin injection characteristics of various 
heterostructures. 

First of all we perform a careful analysis of the regime 
to which GaAs and ZnSe-based structures (most com- 
monly used in spintronics experiments) belong. We point 
out that, with the exception of lightly doped samples 
at room temperature, these systems are generally in a 
regime which cannot be described by non-degenerate 
electron statistics. In particular, the low temperature 
heterostructures used in Refs^andi are in the degener- 
ate regime (see Fig. \Q . We then study these structures 



under the effect of an applied electric field and introduce 
the "semiconductor degenerate regime" characterized by 
the inequalities kgT < £p < e|_B|L s , where Ef is the 
Fermi energy of the carriers in the semiconductor and E 
is the applied electric field. This high-field regime oc- 
curs, for comparable materials and doping densities, at 
fields which are one order of magnitude smaller than in 
the non-degenerate high-field regime. In particular, the 
use of non-degenerate statistics in heterostructures based 
on GaAs, ZnSe or comparable materials, at a few Kelvin, 
can overestimate the effect of the electric field by orders 
of magnitude (see Sec. IV). 

As for Coulomb interactions, we demonstrate that the 
experimentally relevant regime corresponds to the one 
in which the many body effects are at their strongest. 
Neglecting these effects can induce both quantitative as 
well as qualitative errors. Coulomb interactions affect 
the system through the spin Coulomb drag effec1i£*§*£ and 
by reducing the longitudinal spin stiffness. In particu- 
lar we will discuss how the spin drag explicitly modifies 
spin-transport whenever materials with different mag- 
netic properties are joined. 

The numerical part of this paper concentrates on 
sandwiched structures (ferromagnet (FM) /non-magnetic 
semiconductor (NMS)/ferromagnet) which are of the ut- 
most importance in device applications based on the 
magnetoresistance effect. We have performed a detailed 
analysis of the spatial dependence of various quantities 
such as electrochemical potentials, density and current 
polarizations for both parallel and antiparallel alignment 
of the ferromagnets. At variance with the parallel case, 
we find that in the antiparallel case, the magnetoresis- 
tance can be increased by a strong electric field. 

Another important issue we address is the effect of the 
dependence of the spin diffusion constant D s on carrier 

density when /cgT ~ sp (see inset of Fig. E). We will 
show that, due to this dependence, it is possible to greatly 
enhance the local spin current density by carefully choos- 
ing the system parameters. 
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This paper is organized as follows: in Sec. II we discuss 
the system regimes as doping density, temperature, and 
electric field are varied; in Sec. Ill we derive the drift- 
diffusion equation for spin-transport including both the 
Coulomb interaction and the electric field, and discuss 
the behavior of the upstream and downstream diffusion 
lengths; in Sec. IV we discuss the limits of the non- 
degenerate approximation; in Sec. V we consider the spin 
current and the implications of varying the carrier den- 
sity on this quantity; in Sections VI and VII we analyze 
FM/NMS and FM/NMS/FM heterostructures, respec- 
tively; in Sec. VIII we describe the effects of Coulomb 
interactions; finally in Sec. IX we draw our main con- 
clusions. Two appendices, discussing some details of the 
proposed heterostructure-related formalism and propos- 
ing "exact" equations for describing spin transport in 
non-interacting degenerate heterostructures conclude the 
paper. 



II. DOPING DENSITY, TEMPERATURE AND 
ELECTRIC FIELD -RELATED REGIMES 



This work focuses on electrically driven spin transport: 
due to the presence of the external electric field, the 
system is characterized by three different energy scales 
whose relative magnitude determines the system behav- 
ior. Such energies are the Fermi energy £p = h 2 k F /2m 
(which strongly depends on the carrier density), the ther- 
mal energy fc#T and the energy related to the electric 
field e\E\L s , where L s — \JD s t s is the diffusion length, 
D s the diffusion constant and r s the spin-flip time. 

Let us first of all define the regimes the system un- 
dergoes depending on the relation between thermal and 
Fermi energies 

fcgT ^> ep non-degenerate regime (I) 
fcgT <C Ef degenerate regime (2) 
fcgT sa Ep intermediate regime. (3) 

The intermediate regime corresponds to the crossover re- 
gion between the first two, roughly the region between 
the dotted lines in Fig. ^ 

We will concentrate our discussion on NMS in the in- 
termediate and degenerate regimes. Such regimes are of 
utmost importance since they include most of the semi- 
conductor structures used in spintronics related experi- 
ments, such as n-doped GaAs 3 - or Zn .97Beo.o3Sei ones. 
This is illustrated in Fig. ^ where we plot the ratio 
ksT / Ef as a function of the carrier density n for n-doped 
GaAs (upper panel) and (Zn,Be)Se (lower panel). Three 
different values of the temperature are considered (as la- 
belled). The regions corresponding to degenerate and 
non-degenerate systems are specified as well. As can be 
seen, up to room temperatures, for densities as low as 
n ~ 5 • 10 16 cm~ 3 in GaAs orn~ 10 18 citT 3 in (Zn,Be)Se 
the system is at most in the intermediate regime; for 
T = 1.6K, a temperature often used in experiments, the 
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FIG. 1: fcsT/eF vs carrier density n for GaAs (upper 
panel, m = 0.067mo, e = 12) and (Zn, Be)Se (lower panel, 
m = 0.5m , e = 12) parameters, at T = 1.6K, 20K, 300K, 
as labelled. Degenerate and non-degenerate regions are indi- 
cated as well. Inset: D s vs carrier density for T — 1.6K and 
GaAs parameters. 



system remains in the degenerate regime for the whole 
relevant density range. 

Since, as just discussed, in most experimental situa- 
tions ksT ~ Ep, the other relevant information to deter- 
mine the expected system behavior is contained in the 
ratio e\E\L s /£p. 

In Fig. |3 we plot this ratio as a function of density and 
for fields up to 1000 V/cm, for both GaAs and (Zn,Be)Se. 
In both cases we have chosen L s = 2/im. We will consider 
this same value in all the calculations presented 13 . The 
high (e\E\L s ^> Ep) and low (e\E\L s <C £f) field regions 
are explicitly indicated. We notice that, even at very low 
densities, it is necessary to apply fields of the order of 
100 V/cm to enter the high field regime. 

As we will explain in detail, the comparison of the en- 
ergy scales, as done in Fig.^and Fig. [21 is a very efficient 
way to understand in which regime the considered sys- 
tem is and deduce how and if the electric field is going 
to affect the diffusion lengths and the spin-transport re- 
lated quantities. In particular at variance with metals, 
in semiconductors the Fermi energy can be so low that, 
even for moderate electric fields, we can fulfill the con- 
dition e\E\L s 3> Ef ^> ksT. This is what we define as 
"semiconductor degenerate regime" to distinguish it from 
the well-know "metallic" degenerate regime in which 
Ef ~> e\E\L s for any reasonable field. In the "semi- 
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FIG. 2: e\E\Ls / ef vs carrier density n for GaAs (upper 
panel) and (Zn, Be)Se (lower panel) parameters, for E = 
1, 10, fOO, lOOOV/cm as labelled. High and low field regions 
are indicated as well. 



conductor degenerate regime" the electric field heavily 
affects spin-transport even if the system is degenerate. 
In order to illustrate this point, we will analyze in de- 
tail the behavior of the penetration lengths Ld and L u as 
electric field and carrier density change. 

Due to its importance, in this article we will mainly 
discuss n-doped GaAs at T = 1.6K; our results are based 
on the interplay among different energy scales, so that 
they can be readily extended to other semiconductor ma- 
terials. 



III. INTERACTING DRIFT-DIFFUSION 
EQUATION 

First of all we will sketch the derivation of the general 
drift-diffusion equation for spin-transport which includes 
Coulomb interactions among carriers. We refer the reader 
to Ref. for the details. 

We start from the drift diffusion equation for the in- 
teracting spin-a charge current 



3a(r) = X! {^apE + eDapVnp^j 
' ' 



(4) 



where o~ a p is the (non- diagonal) homogeneous conduc- 
tivity matrix of the electron gas and D a p is the {non- 



diagonal) diffusion matrix. Throughout the paper Greek 
indices will correspond to the spin- variables f, J,. We sub- 
stitute Eq. in the generalized continuity equation for 
the spin-density components 



dAn a {r,t) _ An a (f,t) Ara a (r,i) V • j a (r) 



dt 



sf,a 



Tsf,a 



(5) 



where An a (r , t) = n a (r, t) — ria , n a (r, t) is the a density 



component, n« its equilibrium value and T s f t0l is the a 
spin-flip relaxation time. Carefully applying 7 the local 
charge neutrality constraint 



An T (r) = -An;(r) 



(6) 



and considering the regime linear in Vn Q , the general 
interacting drift-diffusion equation is obtained 7 -. 

Finally, imposing the steady state condition 
<9(Artf — Ani)/dt — we obtain the steady-state 
interacting drift-diffusion equation 



(An T - A ni ) 



+H S E ■ V(An T 



hD s V 2 (An T 
AnJ = 0, 



Arn) 



(7) 



where t s — (\/r s t t -\ + l/r s y.|) _1 is the spin relaxation 
time, and the effective interacting mobility and diffusion 
constants are given, for a non-magnetic system^, by 



ero/m 



and 



D, 



Us Q 
— bn 



1 - pu/pd' 



(8) 



(9) 



Here po = m/ne 2 TD is the ordinary Drude resistivity 
m the effective mass of the carriers, S the (interacting) 
static longitudinal spin-stiffness 7 , n = nj the carrier 
density and pn the spin-transresistivity, which measures 
the momentum rate exchanged between spin up and spin 
down carriers^* 7 .. Eq. is valid independently of the 
value of ksT/sF- 

If the sample is homogeneous and E || x the solution 
to Eq. is given by 

An T - An x = Ae x/Lu + Be~ x/Ld (10) 
with A, B determined by the boundary conditions and 



(11) 



L u ^ are the interacting upstream and downstream "pen- 
etration" lengths 12 . Such quantities correspond to the 
average decay lengths of an injected spin-unbalance in 
the upstream (L u ) and downstream (Lj) directions. 

The upper panel of Fig. [3] shows the behavior of 
Ld, as a function of density for GaAs at T=1.6K and 
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FIG. 3: Upper panel: downstream diffusion length Ld vs car- 
rier density for GaAs parameters, T=1.6K and E=25V/cm. 
Solid line refers to the interacting calculation, dashed line to 
the non-interacting approximation and dashed-double dotted 
line to the non-degenerate limit. Inset: Upstream diffusion 
length L u vs density n for GaAs parameters at T — 1.6K and 
three different fields (as labelled in V/cm units). Line types as 
in the main panel. Lower panel: downstream diffusion length 
Ld vs carrier density for GaAs parameters, T=1.6K, 20K and 
300K (as labelled) and E=2500V/cm. Line types as in the 
upper panel. 

E = 25 V/cm. The solid line represents the fully inter- 
acting calculation, the dashed one the non-interacting 
approximation and the dashed double dot line the non- 
degenerate limit obtained by setting D s = [isksT/e in 
Eq. (|lll> . In this, and in all the following calculations, 
we have set fi s — 3 x 10 3 cm 2 /Vsi£ However we stress 
that, at least in the degenerate regime, for which usually 
Pu/pd *C 1, L u> 4 is basically independent from fx a (see 
Eqs. dJ and ©). 

Let us focus on the crossovers undergone by the system 
when the carrier density is increased. At this tempera- 
ture even at very low densities the system is degenerate 
(ef » ksT), so the interesting interplay is between 
the Fermi energy and the electric field related energy: as 
shown in the figure, depending on the system density, 
the Fermi energy can be low enough to enter, even for 
moderate electric fields, the "semiconductor degenerate 
regime" e\E\L s >> £p. The existence of such a regime 
implies that, even for a degenerate system, the drift term 
in Eq. cannot be neglected, i.e. it is not possible to 
describe spin-transport using a diffusion equation for the 



electrochemical potentials, as usually done for metalai*. 
Even in degenerate semiconductor systems in fact, the 
drift term severely modifies the penetration lengths L Uf d, 
varying their values over order of magnitudes (as clearly 
shown for Ld in Fig.[3J). As a rule of thumb, we see that in 
this regime the order of magnitude of such a variation is 
set by the ratio e | E | L s /ef- The "semiconductor degener- 
ate" regime extends to higher densities when the applied 
field is increased. Notice that, even in this regime, the ac- 
tual value of Ld strongly differs from its non-degenerate 
limit. By increasing the carrier density, the system will 
eventually enter the "metallic" regime (as indicated in 
the figure), in which the drift term can be neglected and 
L u ,d recover their unperturbed value L s . 

Similar results, i.e. order of magnitude variations de- 
pending on the carrier density and quantitative impor- 
tance of Coulomb interactions are obtained for the up- 
stream penetration length L u (see upper panel inset), 
with the key difference that in the "semiconductor de- 
generate regime" L u is strongly reduced by the electric 
field effect, being the product L u ■ Ld — L s constant. 

It is interesting to ask what happens in the high field 
regime e|i?|-L s >> [ep, ksT) when both temperature 
and carrier density varies. In Fig. [3 lower panel, we plot 
Ld in the high field regime (E=2500 V/cm), as a function 
of carrier density and for three different temperatures, 
T = I.6K, T = 20K and T = 300K. In addition to the 
crossovers discussed in relation to the upper panel, for 
low densities and T — 300K the system enters the non- 
degenerate regime ef « ksT (see also Fig. HJ. In this 
region in fact, the non-degenerate approximation coin- 
cides with the interacting calculation. Fig. [3] shows that 
the greatest variations due to the presence of the electric 
field are actually reached by lowering the temperature 
and for any reasonable density and the materials con- 
sidered, correspond to the "semiconductor degenerate" 
regime. 



IV. LIMITS OF THE NON-DEGENERATE 
APPROXIMATION 

As demonstrated in Sec. II, the use of non- 
degenerate electron statistics to describe carrier dynam- 
ics in commonly-used semiconductors is justified only for 
very lightly doped semiconductors at high temperatures. 
For this reason, when analyzing spin transport, particu- 
lar care must be taken in determining the correct system 
regime: for some other problems an easy-to-implement 
approximation, used outside its supposed validity regime, 
can still yield reasonable estimates^; unfortunately this 
is not the case for the non-degenerate approximation in 
semiconductor spin-transport. It is tempting, since it 
leads to simpler analytical formulas, to approximate D s 
in the various quantities (from the penetration lengths 
to the spin current, from the injected density polariza- 
tion to the magnetoresistance) , by its non-degenerate, 
classical expression. Unfortunately as it has been ex- 
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perimentally shown£ and theoretically confirmedii, the 
spin diffusion constant varies over orders of magnitude 
going from the non-degenerate to the intermediate and 
degenerate regime (see inset of Fig. ^ hi which D s is 
plotted as a function of n for GaAs parameters). This 
implies that approximating spin transport behavior by 
non-degenerate expressions in the wrong regime can lead 
to quite large errors. In Fig. |31 Fig. 0] and Fig. ^3 wc 
have compared the interacting calculation (solid line) 
with its non-degenerate approximation (dashed double- 
dot line) for some of the most relevant quantities, as pen- 
etration lengths, spin current, and magnetoresistance re- 
spectively. At low temperatures, even for densities as low 
as ij ~ 5 x 10 16 cm -3 and fields as small as E=10V/cm 
(low-field regime), this approximation can overestimate 
the field effect by orders of magnitude. Such a discrep- 
ancy is quite general and is found also for quantities such 
as injected density and current polarizations (not shown). 



V. THE SPIN CURRENT 

We will now focus on the spin-current in a non- 
magnetic system. Starting from the drift-diffusion 
Eq. (0J, and imposing the conditions D s = _Djj — = 

D U - D lt and = A*TT + MU = + MJT which are 
obtained from the general formula for D s and fj,^ in the 
non-magnetic limit and in the regime linear in Vn Q , the 
charge current becomes 

j{r) := j T (f) + j x (r) = eEfi s ( n] + n L ) (12) 

while the spin current assumes the form: 

l(r) := j T (r) - j^r) = jP(x) + eD s nVP{x) (13) 

with P(x) = (n-f — ni)/n the density polarization. If we 
consider a spin unbalance injected at x = 0, i.e. such 
that in Eq. fijfy either A = (x > 0) or B = (x < 0), 
An-f — An^ will acquire a simple exponential form; the 
upstream (x < 0) (downstream (x > 0)) spin current will 
then be given by 



J 



' u ' d {r) 



J T (-e) 7 — x 



P(x) = (j+?/)P(x). (14) 



Eq. i|14fl is very interesting because it clearly shows how 
the spin current is composed: results in fact as the 
sum of two distinct parts, the first corresponding to the 
total drift current j and the second being a total diffu- 
sion current jo- Both the components are "weighted" by 
P(x), the percentage of spin-polarized carriers. It is in- 
teresting to notice that in the downstream case, j \\ jo, so 
that both diffusion and drift phenomena positively con- 
tribute to the spin-current, while in the upstream case 
j and jf) are in competition. It is easy to see that in 
the limit of very large electric fields, 1/Ld — > so that 
Jd (f) = jp( x ), while l/L u -> \E\fi s /D s and jf{r) ex- 
actly vanishes. Keeping in mind that the diffusion con- 
stant increases by orders of magnitude when going from 
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FIG. 4: Diffusion current density jr> and total current density 
j (rescaled by total carrier density n) vs electric field E in the 
downstream direction for GaAs parameters, T=1.6K and two 
different carrier densities (as labelled). Solid lines refer to 
the interacting jo/n, and the dashed-double-dot line to the 
non-degenerate limit of Jd/ti. Dotted line refers to j/n. Inset: 
correction to the non- interacting approximation (j l D — j'S*)/iS' 
vs carrier density for E — 2.5, 25, 250V/cm. 



the non-degenerate to the degenerate regime, Eq. lfTl|) 
suggests that a method to increment the spin current for 
a given electric field should be to enter the degenerate 
regime (for example by lowering the temperature), in- 
creasing in this way the diffusive part of the current. If 
instead the doping (or the temperature) is kept constant, 
an increase in the electric field will increase the drift cur- 
rent j. In this last case though, 1/Ld will decrease, and 
to the incremented j will correspond a decreasing of j^. 

We illustrate this behavior in Fig. 0] In this figure 
we plot jf) and j at T = 1.6K rescaled by the trivial 
(and common) dependence on the carrier density n as 
a function of the applied electric field and for two dif- 
ferent densities (as indicated), j/n (dotted curve) does 
not depend on the density nor on Coulomb interaction. 
The fully interacting jjj/n corresponds to the solid line. 
If we fix the density and increase the electric field, jo 
decreases and j increases; but if we fix the electric field 
(considering a moderate electric field), and increase the 
density, jr> increases while j does not change, so that 
there is no competition between the two quantities. We 
underline once more that this is not due to a trivial lin- 
ear scaling with the carrier density, but to the behavior 
of the diffusion constant when entering the degenerate 
regime. Regarding P(x), the second factor in Eq. (|14fl . if 
we consider as injector a not-fully polarized metal in the 
degenerate regime and for moderate fields, we can easily 
deduce from Eq. (|32|l (next session), that P(0) ~ n 1 / 3 , 
i.e. P(x) slightly increases with carrier density. This im- 
plies that in such a regime our conclusions hold for the 
full spin current j*. 
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FIG. 5: Schematic view of the FM/NMS (left) and 
FM/NMS/FM (right) heterostructures considered. 



VI. FM/NMS JUNCTION 

We will now concentrate on FM/NMS heterostruc- 
tures. The schematic view of the system is illustrated 
in Fig. [5^. We will first derive explicit expressions for 
the quantities of interest (electrochemical potential, den- 
sity and current polarization, magnetoresistance, etc) in 
the interacting case, for general carrier density n and 
temperature T and then concentrate on the degenerate 
regime. As underlined before, this regime is experimen- 
tally relevant, and, as discussed in detail in Sec. VII, in 
it Coulomb interactions are quantitatively important. 

In the non-interacting approximation it is possible to 
write explicit expressions for the chemical potential in 
both the non-degenerate and degenerate regime. By us- 
ing them, the equations for the quantities of interest can 
be derived (see respectively Ref. and Appendix B). 
Since we want to derive instead equations which include 
Coulomb interactions and are valid independently from 
the doping regime, let us consider the Taylor expansion 
of the interacting chemical potential 



chem/ \ chem , 

M<7 { n h n l)-^Qa + 



chem , 9Mf eW A „ , Ql4 



chem 



-An n 



dn s 



-An* 



If we use the relation 



Sa6 



chem 



dm 



(15) 



(16) 



between the {aS} element of the spin-stiffness tensor and 
the chemical potential, and the local charge neutrality 
constraint Eq. JBJ, we can rewrite such a potential as 



chem chern 
AV ~ H'Orr 



{Saa S aa )An 



(17) 



The limits of this approximation are discussed in Ap- 
pendix A. The excess electrochemical potential associ- 
ated to Eq. l(T7|) then becomes 



« (So 



chem 



S aS )An a + eE-x + B. 



(18) 
(19) 



Eq. (|19f) is general and is valid both for polarized and 
non-magnetic materials. 



By using the relation between the spin stiffness and 
the spin stiffness tensor components 



S - (Stj - S n + Sn - Sit)/ 4 , 



(20) 



the symmetries in the NMS (S^ — S-\i = Sn — S^), and 
the local charge neutrality constraint, the electrochemical 
potential in the NMS (x > 0) is given by 

A/i T(i) = ±SnP(x) + eE -x + B 

= ±^(l-p n /p D )P(x)+eE-x + B{21) 

Ps 

where the plus (minus) sign refers to the | (J.) spin com- 
ponent and we have used Eq. @ in the second expression. 

If we assume that the FM (diluted magnetic semicon- 
ductor or metal) is in the "metallic" degenerate regime 
(as it is often the case), the electrochemical potential on 
the FM side {x < 0) will be2J& 



( A/i T 
1% 



J ■ x 



Cej • x 




In this expression Coulomb interactions are included 



through o~l = al a + 0^5™- Since in this regime the 
spin transresistivity is much smaller than the Drude 
resistivity^, we will neglect the former in our calculations. 

If we assume transparent interfaces, the boundary con- 
ditions to be satisfied at x = will be 



3s\q- — 3s\o+ 

a ^tI - = a mtI + 

A Mll - = A AHl +- 
The spin current components are given by 



EdAfis 
O-aS— X, 

5 



Ox 



(23) 
(24) 
(25) 



(26) 



and the spin currents in the two materials will then be- 
come 



3 



2C , x , 

- MP/ + Z7 exp( Z7- 



jr = p(x)u+j d D ) 



(27) 
(28) 



where the indexes /, sc correspond respectively to the 
FM and to the NMS, p f = (ct( - af)/(a^ + of) and % 
is defined by Eq. lfTI)l. 

By using Eqs. (|21|I - (|28JI . after some straightforward al- 
gebra, we derive the expressions for the quantities of in- 
terest. 

The electrochemical potential in the NMS will be 
Afi m) = -SnP(0)\p f T exp(-x/L d )} +eE-x, (29) 

where the - (+) sign refers to the | (J,) spin component. In 
particular the expression for the interface band splitting 
is given by 



A/x T (0)-A M; (0) = 25nP(0) 



2Snp f (30) 
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We see that the splitting saturates for high electric fields. 
The relation Sn = D s e{\ — Ph/pd)/p>s implies on the 
other hand that this limit can be increased by increasing 
the diffusion constant D s (for example by increasing the 
system density as previously discussed). 

The magnetoresistance (which is related to the elec- 
trochemical potential sum at the interface and will be 
discussed in detail later) can be written as 



A MT (0) + A Mi (0) _ SnP(0)p f 



2e|j| e\j\ 
The equation for the injected density polarization is 



(31) 



Pf 



""" =L - { i-hj 



while the one for the injected current polarization is 

Pf 



a(0) 



L u(l - pu/pd) + i 

(32) 
(33) 



^L u (l-p n / PD ) + l 



where we have used the relations Eu s /D s — (1/L U — 
1/L d ). We stress that Eqs. lf25 |) -l|53 )) are valid in all den- 
sity and field regimes and represent the fully interacting 
expressions. 

Fig. illustrates the behavior of the j and J, elec- 
trochemical potentials (excluding the trivial linear de- 
pendence on i, as we will do in all figures related to 
such quantities) as a function of x across the FM/NMS 
junction and for increasing field and FM polarization. 
The parameters used in the calculations we present are 
n = 5.2- 10 16 cm" 3 (e F = 7.6meV), a f /a sc = 95, 
Lf = 20nm, L s = 2/zmi£. In the upper panel of Fig. 
the FM polarization is pf = 0.5. As the figure shows, the 
splitting of the FM sub-bands at the interface increases 
asymmetrically in respect to their asymptotic values (at 
— oo), with the electric field due to the difference be- 
tween the conductivities <r^ , erj . In the inset we plot 
the same electrochemical potentials but on the NMS side 
of the junction (x > 0). As expected the potentials are 
now symmetric with respect to their asymptotic value 
(at oo) and the extension of the spin-polarized region 
in which V/i CT ^ const extends in space as the electric 
field increases. As can be seen, the potential drop at the 
interface (the difference between the values of the elec- 
trochemical potentials at ±oo), which is related to the 
magnetoresistance, increases with the field. 

The lower panel of Fig. illustrates the influence of 
increasing the FM polarization on the electrochemical 
potentials at the junction and in the high field regime 
(E = lOOOV/cm). Fig.[7|shows in detail the behavior of 
the related quantities as pf is increased from 0.5 to full 
polarization. While most of the considered quantities 
(from the injected current and density polarizations to 
the sub-band splitting to the magnetoresistance) mono- 
tonically increase with the FM polarization to saturate 
when pf = 1, the | sub-band split (dashed line labelled 
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FIG. 6: Upper panel:Up and down electrochemical poten- 
tials (without the term linear in x) at the junction vs x for 
three different fields (as labelled) and pj = 0.5. Inset: as 
in main panel but on the NMS side of the junction. The 
non-interacting approximation (dashed line) is plotted as well. 
Lower panel: As in the upper panel but for three different FM 
polarizations (as labelled) and E — lOOOV/cm. 



with |) increases at first to rapidly drop to zero when the 
FM polarization saturates. 

As in non-degenerate systems^, Pq saturates towards 
the ao value as the field increases (not shown). 



A. Magnetoresistance 

A very important quantity for building "on-off" type 
devices as magnetic field sensors or read heads for mag- 
netic hard disk drives^ is the so called magnetoresistance 
R m . It can be defined as the part of the equivalent re- 
sistance R eq of the device which depends on its magnetic 
properties. 

For the considered structure R 



R, 



_ A^(A SC ) 



eq is given by 
Am(-A/) 



eq 



e\j\ 



(34) 



where ^ sc (f) ar e the device lengths in the NMS and FM 
regions, and A/i = (A/ij + A/i|)/2 indicates the variation 
of the electrochemical potential with respect to the situ- 
ation in which no electric field is applied, and no junction 
is present. We will define R m as R eq minus the trivial, 
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FIG. 7: Up and down electrochemical potentials (dashed 
lines), their difference (dotted line), the current and density 
polarizations (solid line) calculated at x — as well as mag- 
netoresistance R m (dashed-dotted line) vs FM polarizations 
p f for E = lOOOV/cm. 



non-magnetic device resistance Rt r when A sc 3> L s and 
\f > LfUL For a FM/NMS junction we obtain 



Rr, 



L, 



L, 



a sc (i s \E\t, 



(1 



PU 

PD 



)p f P(0). 



(35) 



Eq. I|35(l is very interesting: it clearly shows that R m de- 
pends on the intrinsic resistance of the NMS L s /a sc , on 
the ratio between diffusion and drift lengths L s / h s \E\t s 
and it emphasizes the dependence of R m on the spin drag 
effect. It also shows that the magnetoresistance is pro- 
portional to the probability p/P(0) that a spin crosses 
the interface without loosing its spin polarization, i.e. it 
preserves its magnetic properties when entering the NMS. 
Eq. i|35|) shows that, whichever density regime is consid- 
ered, the electric field destroys the magnetoresistance. 

In Fig. ( lower panel inset, curves labelled "ljct") 
we plot the ratio between R m and the total device 
resistance R eq in respect to the applied field and for 
^f(sc) — 10^/(sc) : such a ratio increases by orders of 
magnitude with the FM polarization pf. At £'=0.01 
V/cm, R m = 0.28n£!cm 2 forp/ = 0.5 and R m = 8fiilcm 2 
for pf = 1. 



VII. FM/NMS/FM STRUCTURE 

We will now consider the tri-layer structure of Fig. [5J3. 
According to Eq. (|10|l and to the local charge neutrality 
constraint Eq. (JBJ, in the NMS the excess spin density 
components are given by 



A? 



'T(l) 



= ± 



a r x n a 1 ( x - x o) X 

A exp (- — ) + A x exp ( ^— '- ) 



(36) 



where Xq > corresponds to the second interface and 
the + (-) sign to the f (J,) component. The chemical 
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FIG. 8: Upper panel: Up (solid line) and down (dashed line) 
electrochemical potentials (without the term linear in x) vs 
x for three different fields (as labelled) and pf = 0.5. The 
two ferromagnet polarizations are parallel. Inset: as in main 
panel but for E= 10000 V/cm. Lower panel: As in the upper 
panel but with antiparallel ferromagnet polarizations. 



potentials in the three materials are now 
j ■ x l 1 



A MT(i) 



1 



+ Cz,ej ■ x 




ej ■ x 



+ C R ej ■ x 



where the indices L, R stand for left and right ferromag- 
nets and the + (-) sign in Eq. (55)1 corresponds to the | (J.) 
component. The six constants {Aq, A\, B sc , Bf, Cl, Cr} 
will be determined by boundary conditions similar to the 
ones for the single interface, i.e. 



A ^Tlo-,, 



a mtI +,, 

A Mllo+, 3 



(40) 

(41) 
(42) 



First of all let us consider the spin density and current 
polarizations. In order to write the expressions for these 
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quantities we need to derive the expressions for Aq, A\. 
For general polarizations p^ L , p R , after some algebra, the 
expressions for Aq, A\ can be written as 



i 



l\ — CXp( — Xq 



l R 
"1 -1 



(43) 



M = -2 Ld lr u 



where 



1 



d%-exp(-x [— + — ))d L ^ 



L R,L 



= 1±G 

= l±G RtL 



L d 

R„L — T 



R,L 
J R,L 



(44) 

(45) 
(46) 
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Grx = 



^[1-(P^) 2 ](1-^)- 



(47) 



The spin density and current polarizations are then ex- 
pressed as 



P(x) 
a{x) 



A exp(- — ) + A x exp( v " ) 

J-'d J-'u 



2 
n 

2/1 1 

n \L U Ld 

, X . 1 Ax — Xq) . 1 

A exp(- — )— - A\ exp( l o; )-H 

J^d Li u Li u L d 



(48) 



49) 



In the interesting case in which p^ L = ±p R — p^ (parallel 
and anti-parallel polarization of the two FM's, which are 
considered, for the rest, equivalent) the current polariza- 
tion at the two interfaces becomes: 



a(0) = 



P 



d+u+ - exp(-a;o + ±) )d L u] 



d R ~ exp(-xo I — + — I )u R 



1 

L d 



a(x ) 



±exp(-^-)^(L„ + L d ) 
±pf 



d+u+ - cxp(-x + -fc))d L u R 



(50) 



1 



1 



exp(-a; ( — + — ))d 



, / x \ GR it , r \ 

±exp(- — )— T [L U + L d ) 
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(51) 



FIG. 9: Upper panel: density (solid line) and current (dot 
line) polarizations vs x for E=10000 V/cm and pj = 0.5. 
The two ferromagnet polarizations are parallel. Inset: as in 
main panel but for three different fields (as labelled in units 
of V/cm). Lower panel: density (solid line) and current (dot 
line) polarizations vs x for E=10000 V/cm and Pf = 0.5. The 
two ferromagnet polarizations are antiparallel. Upper inset: 
as in main panel but for E=l V/cm. Lower inset: as in main 
panel but for E=100 V/cm. 



where the upper (lower) sign refers to the parallel (an- 
tiparallel) configuration. 

Finally the expressions for the other constants which 
appear in the electrochemical potentials are given by 



B sc 
Bf 

C L 

Cr 



-SnP(Q)p f L 



(52) 



f 

R 



a 



\E\ 

lil 



Sn 



xo + ^\p f L P(0)-p R P(x )](53) 



2 e\j\ 
Sn P(xq) 



Pl ) 



'lil 



4(i -p f R 2 )- 



(54) 
(55) 



In Fig. |H| we plot the behavior of the f (solid lines) 
and i (dashed lines) electrochemical potentials in the tri- 
layer structure for the interesting case in which the NMS 
width xq is much smaller than the diffusion length (xq — 
L s /10). The upper panel shows the electrochemical po- 
tentials when a parallel configuration {pjL = PfR = 0.5) 
is considered. The lower panel shows the two junction 
system when the polarizations of the two FM's are op- 
posite (anti-parallel case,p/L = —PfR. = 0.5 ). For both 
configurations, as the field is increased the potential drop 
due to the presence of the NMS increases. 



10 



We notice immediately that in the parallel case the 
two spin components of the electrochemical potentials 
cross. This crossing was predicted to occur in Ref. [l| 
and is necessary in order to obtain a different potential 
drop for the | and | spin components and hence a finite 
current polarization even at low electric fields. A clear 
effect of the electric field enhancement is to "push" the 
crossing between "f and | components toward the second 
interface by increasing the penetration length Ld rela- 
tive to the first interface and decreasing the upstream 
penetration length L u relative to the second interface. 
This behavior becomes extreme for very high fields (see 
inset where E = lOOOOV/cm), where the electrochemi- 
cal potential components are almost parallel until they 
suddenly switch before the second interface. 

In the antiparallel configuration (lower panel), while 
the potential drop across the whole structure behaves 
similarly to the parallel case, the electrochemical poten- 
tials do not cross and this behavior destroys the spin cur- 
rent at very low fields 1 . The electrochemical potentials 
present instead a quite distinct feature: for increasing 
fields a spin accumulation builds up for iwioi just be- 
fore the second interface, underlining the increasing resis- 
tance that spins encounter in crossing this interface when 
the field (and with it the polarization inside the NMS) 
increases. This anticipates that in the antiparallel con- 
figuration the magnetoresistance is expected to increase 
with increasing applied electric field. 

Fig. presents the behavior of the density (solid lines) 
and current (dotted line) polarizations as a function of x 
for both the parallel (upper panel) and the anti-parallel 
(lower panel) configuration, when the field is increased 
from the low field regime (E = lV/cm, E = lOV/cm), 
to the intermediate (E = lOOV/cm), to the high field 
one (E = lOOOV/cm, E = lOOOOV/cm). The main panel 
presents the highest field situation. In both configura- 
tions, as expected from the behavior in the single junc- 
tion system, for low to intermediate fields -P(O) < ct(0), 
while -P(O) — * a(0) in the high field regime. On the other 
side, in the parallel configuration, while, due to the choice 
xq <C L s , a{x) is basically constant even at very low 
fields, the density polarization P{x) switches sign when 
the two electrochemical potentials cross. This inversion 
in the polarization of the spin population unbalance is 
pushed toward the second interface as the field increases. 

The situation is different in the antiparallel case in 
which, at low to intermediate fields, an almost constant 
P{x) corresponds to a constant drop of the current po- 
larization a(x) as a function of x. For very low fields 
(E = lV/cm), and just before the II junction, it be- 
comes even slightly negative, as could be expected by 
simply considering the diffusion current at this interface. 
The drop in a(x) is due to the difficulty in establish- 
ing a spin current through the system due to the abrupt 
decrease of the |-spin density of states in the second fer- 
romagnet. When the field is increased, the decrement of 
the current density becomes less drastic while a spin ac- 
cumulation starts to build up at the second interface (see 



lower panel of the inset). In the high field regime finally, 
it is possible to establish a basically constant spin density 
current through the device, at the price of building up a 
strong spin density accumulation just before the second 
interface (main lower panel). 



A. Magnetoresistance 

Let us now concentrate on the magnetoresistance. R eq 
in the FM/NMS/FM structure is given by 

R eq . M^. + A/.«)-M-A/,l) | (56) 



e\j\ 



where the NMS extends from to xq. By writing 
^f,R(L) — n LfR(L) and using for the electrochemical po- 
tentials Eqs. lj37|) . il3S|l and we obtain 



Re 



L 



fR 



OfR 



LfL 

<JfL 



X 



~ exp(-n)) ' 

O sc Vs\E\T s 
■\p f LP(0)-pf R P(x O )} 



( l_£Ti) 

PD 



(57) 



The first line is clearly the trivial, non-magnetic term. 
For n> 1, R m is given by 



Rm. — 



<T SC fJ, s \E\T s p D 

■\p fL P(0) - PfR P(x )} 



(58) 



The above equation shows that R m (i) is proportional 
to the intrinsic resistance of the paramagnet L s /a sc ; (ii) 
it depends directly on the ratio between the diffusion 
length L s = \/D s t s and the drift length fj, s \E\T s ,i.e. it 
tends to diminish with increasing field since the spins 
will be drifted longer and less resistance occurs (iii) it 
depends directly on the spin drag which enhances it; (iv) 
it is proportional to the difference p/z,P(0) —pfnP(xo) = 
\pfLP(0)\ + \pfnP{xo)\. This term can be seen as 
the probability that magnetic properties remain relevant 
through the whole structure. It implies that, as long as 
the products \p/P\ are not vanishing, magnetoresistance 
is accumulated at the different junctions . (v)In the im- 
portant case in which p{ = ±p R = p? , the last factor 
becomes pf[P(0) =F P(%o)]'- this shows that in the anti- 
parallel case the magnetoresistance is proportional to the 
sum of the density polarizations just after the first inter- 
face and before the second one. This implies that it is 
enhanced by the spin accumulation at the second inter- 
face due to the increasing of the electric field (see Fig.|5J|. 

Fig. ^]] illustrates the behavior of the magnetoresis- 
tance R m as a function of the electric field. The main up- 
per panel refers to the parallel configuration. We clearly 
see that the field decreases the magnetoresistance. The 
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FIG. 10: Upper panel: magnetoresistance R m for the parallel 
configuration vs electric field for Pf = 0.5. Its non-degenerate 
approximation (dashed-double-dot line) is plotted for compar- 
ison. Inset: as for main panel, but for the antiparallel config- 
uration. Lower panel: magnetoresistance variation from the 
the parallel to the antiparallel configuration vs electric field 
for pf = 0.5. Its non-degenerate (dashed-double-dot line) 
and non-interacting (dashed line) approximations are plotted 
for comparison. Inset: ratio Rm/Req vs E for single (ljct) 
and double junction (2jct, parallel configuration) structures 
and different FM polarizations (as labelled). For pf — 1 the 
dashed-double dot line indicates the non-degenerate approxi- 
mation. 



inset shows the corresponding situation for the antipar- 
allel case: as anticipated, in this case the magnetoresis- 
tance is increased by the applied electric field. The lower 
panel presents the ratio of the difference between the par- 
allel and the antiparallel magnetoresistances to 
Rfn- This relative change in magnetoresistance, respon- 
sible for the functioning of spin-valve based devicesii, 
is increased by the field effect. The non-degenerate and 
non-interacting approximations are shown for compar- 
ison. The inset presents the ratio R\ n lR v eq (curve la- 
belled "2jct") as a function of the applied field. It is 
interesting to notice that this ratio is about two or- 
ders of magnitude higher than the single junction result 
with the same FM polarization. In systems such that 
L sc /o~ sc 3> Lf/o~f in fact, while the magnetoresistance 
R m per junction does not change significantly between 
the single and double junction systems, the non-magnetic 
component of R eq is mainly determined by the NMS layer 
length, which strongly decreases in the the tri-layer de- 
vice. This suggests that, an important requirement for 
optimizing magnetoresistance-based devices, would be to 
minimize the NMS layer length. 



VIII. COULOMB INTERACTION EFFECTS 

Coulomb interactions enter the drift-diffusion equation 
(Eq. 10) through the diffusion and mobility constants; in 
the non-magnetic limit (Eqs. (JSJ and JDJ), through the 
diffusion constant D s only, by affecting the spin-stiffness 
<S£ and generating the spin transresistivity p^j^Si 

First of all, we want to point out that the experimen- 
tally important regime ksT ~ ep (see Fig. 0, is the 
same in which Coulomb interactions between carriers of 
opposite spin become relevant. In fact for ksT ~ ep the 
spin transresistivity reaches its maximum and can be- 
come of the same order of the Drude resistivity^, while for 
ksT ~ ep the spin-stiffness displays the maximum devia- 
tion from its non-interacting approximation™. These two 
combined effects reduce the non-interacting approxima- 
tion to D s even by 50%, depending on carrier density and 
temperaturat According to the functional dependence 
on D s of the various quantities of interest, Coulomb in- 
teraction effects can get amplified, as we will explicitly 
show for the downstream spin penetration length and for 
the spin current. 

The strong contribution of Coulomb interactions to 
Lu.d values is underlined in the upper panel of Fig. 1111 
(as well as by the direct comparison between solid and 
dashed lines in Fig. in which we plot, for GaAs at low 
temperature, the correction {L l d — L™)/L™ as a function 
of carrier density and for different fields. Here and in the 
following i,ni stand for interacting and non-interacting 
respectively. As the figure shows, the correction can 
be of the order of 100%, definitely not-negligible when 
quantitative calculations are required. This correction 
becomes the largest in the "semiconductor degenerate 
regime" 15 e\E\L s » ep » ksT : when the penetration 
lengths approach their high field limits Ld ~ L? s E[i a j 'D s , 
L u sa D s /(En s ), which strongly depend over the diffu- 
sion constant. At low densities such a limiting behavior is 
reached already for E = 25V/cm, while the curve corre- 
sponding to E = 2.5V/cm shows smaller corrections. For 
higher densities the limiting behavior is reached at larger 
fields (see the curve corresponding to E — 250V/cm), 
though the interacting correction remains significative 
even for fields as low as E = 25V/cm. 

The lower panel of Fig. ^] shows the importance of 
Coulomb corrections in the high field regime, when tem- 
perature and carrier density are varied. The largest cor- 
rections are present at T = 1.6K and low densities and 
are due to the spin stiffness reduction. Even at temper- 
atures as high as room temperature though, Coulomb 
interactions remain relevant and the correction can be 
as high as 30%. We underline that in this last case, it 
is mainly due to the spin Coulomb drag (SCD) effect 5 ^ 
which enters the diffusion constant through the spin tran- 
sresistivity (see Eq. (|!JJ|). We want to stress that Coulomb 
interactions modify the behavior of L u c i even qualita- 
tively: as shown for Ld in the lower panel of Fig. |3| 
(T = 20, 300K), due to Coulomb corrections, the non- 
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FIG. 11: Upper panel: \L d — L™\/L™ vs carrier density for 
GaAs parameters, T=1.6K and three different electric fields 
(as labelled). Inset: AR m /R m = (R l m - R™)/R™ vs electric 
field forp/ = 0.5 andp/ = 1 (as labelled). Lower panel: \L d — 
U d \jI7 d vs carrier density for GaAs parameters, T=1.6K , 
20K and 300K (as labelled) and E=2500V/cm. 



degenerate limit is approached from above instead than 
from below when n is decreased. This implies that, even 
if for high enough temperatures the system can actu- 
ally enter the non-degenerate regime, the maximum of 
Ld (the minimum of L u ) is in any case reached outside 
such a regime. 

The inset of Fig. 0] shows the Coulomb correction to 
jf) as a function of carrier density and for three different 
electric fields (as labelled in units of V/cm). Again we 
notice that the functional dependence of jf) on D s am- 
plifies interaction effects which, depending on the carrier 
density, are quite substantial. 

Let us now consider how Coulomb interactions should 
affect heterostructures between materials with different 
magnetic properties. Apart from entering the equations 
of motion through quantities as the diffusion constant or 
the penetration lengths, we expect the SCD to play a 
special role and to affect independently all the quantities 
related to interface properties. This is due to the very 
nature of the SCD effect: in fact it tends to relax spin 
current opposing relative motion between up and down 
spin components, which include opposing all their rel- 
ative variations connected to injection of a drifting spin 
unbalance. This picture is confirmed by all our equations 
related to heterostructures, starting from the expression 
for the electrochemical potential Eq. 1)2 l[l. in which the 
spin Coulomb drag enters explicitly through the factor 
(1 — p^/pu). As a consequence, all quantities derived 



from this potential as injected spin density, current den- 
sity, and magnetoresistance, will explicitly display such 
a factor. 

The spin-transresistivity p^ is a negative quantity, 
so the factor (1 — Ph/pd) will oppose spin injection in 
Eqs. (j32l and while it will increase the magnetore- 
sistance in Eqs. 1)31(1 and l|58|l. Since for some systems 
Pll ~ pe£, the SCD could actually strongly affect the 
magnetoresistance. In the calculations here presented 
though, we have considered only degenerate systems, in 
which the SCD is usually negligible. 

In the inset of Fig. H3 for E — lOOOV/cm, we plot 
both the interacting (solid line) and the non-interacting 
(dashed line) calculations of the electrochemical poten- 
tials: as early discussed, at high fields, interactions in- 
crease the spin penetration length Ld in a noticeable way; 
this correspondingly increases the extension of the polar- 
ized region in the NMS. 

For all calculations related to heterostructures, we have 
considered the carrier density n — 5.2 x 10 16 cm~ 3 : as 
can be expected by the upper panel of Fig. El the in- 
fluence of Coulomb interactions should then be not neg- 
ligible and of the order of 20-25%. This is confirmed 
by the detailed analysis of the various quantities of in- 
terest: in a FM/NMS heterostructure the influence of 
Coulomb interactions on a(0) and P(0) is of the order 
of 25% (see inset of Fig. EJ|. For E -> 0, the ratio 
(P(5 - Pq 1 )/Po 1 -> (Pi - Df)/Df, while as the field is 
increased well inside the high field regime, L u strongly 
decreases and the Coulomb related term in Eqs. (|32|l and 
1)33)1 becomes less important, so that the Coulomb cor- 
rection starts to decrease. For the same structure, the 
inset of Fig. ^2 presents the Coulomb interaction effects 
on R m for increasing electric field and polarization. Here 

ARm/Rm = (Rtn ~ R m)/ R m- We see that such effects 
become important as both field and pf increases. 

When we consider a double-junction heterostructure, 
the effect of Coulomb interactions remains sizeable. In 
Fig. El we plot the interacting (solid line) and non- 
interacting (dashed line) spin density polarization P(x) 
across the NMS, both in the parallel (labelled by "p") 
and antiparallel (labelled by "a") configuration for E = 
lOOOV/cm. Again, as for the one-junction case, the in- 
teraction correction is of the order of 20-25%. Similarly 
from the lower panel of Fig. El we see that Coulomb in- 
teractions at high fields and at the considered density can 
influence the ratio (Rfn — R^/R^ by the same amount. 



IX. CONCLUSIONS 

We have presented a detailed study of spin injec- 
tion into degenerate semiconductors, analyzing the com- 
bined effects of applied electric field, carrier density and 
Coulomb interactions. 

We have clarified the interplay of the different en- 
ergy scales involved in the problem, proposing an effi- 
cient way to understand the regime of the system and 
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FIG. 12: Double junction system: Density polarization P(x) 
vs x for the parallel (labelled as "p") and antiparallel ("a") 
configurations (solid line). The non-interacting approxima- 
tion (dashed line) is plotted as well. Inset: Single junction 
system: AP = (Pi - PS*) /PS 4 and Aa„ = (aj - af)/af 
vs electric field. 



estimate the effect of the applied electric field. In par- 
ticular we have focused on ferromagnet/non-magnetic 
semiconductor and ferromagnet/non-magnetic semicon- 
ductor/ferromagnet structures, carefully comparing the 
behavior of the relevant quantities (electrochemical po- 
tentials, spin density, magnetoresistance...) when the po- 
larization of the second ferromagnet is switched from par- 
allel to antiparallel or its absolute value and the electric 
field are increased. 

We have pointed out, in the anti-parallel case, that the 
presence of a spin-current through the device is related to 
the onset of a spin accumulation at the second interface, 
which in turn increases the device magnetoresistance: at 
variance with the parallel case, we demonstrate that in 
the anti-parallel configuration the magnetoresistance can 
be increased by a strong applied electric field. Addition- 
ally we show that this field increases as well the relative 
change in magnetoresistance |i?f„ — Rf^/R^, exploited 
in spin-valve based metal devices^. 

We have proposed a general formalism which includes 
Coulomb interactions and applies to FM/NMS junctions 
independently of their carrier density regime. In partic- 
ular it allows also to study heterostructures in the in- 
termediate regime e ss fcgT. "Exact" expressions for the 
case in which a non-interacting degenerate system is con- 
sidered have also been provided (see Appendix B). 

Our calculations show that, in the experimentally im- 
portant regime ej ~ fcsT, corrections due to Coulomb 
interactions are relevant. They also demonstrate the im- 
portance of carefully identifying the system regime (for 
example by comparing energy scales as we suggest in 
Sec. II): non-degenerate approximations in fact, when ex- 
tended to the wrong regime, can lead to errors of orders 
of magnitude. 
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APPENDIX A: LIMIT OF THE CHEMICAL 
POTENTIAL APPROXIMATION 

The approximation we proposed for the chemical po- 
tential, Eq. (|17fl . can be thought at first sight to be valid 
only in the low polarization limit P(x) -C 1. 

In this appendix we want to estimate the limits of this 
approximation in the important case of a degenerate sys- 
tem. 

In a non-interacting system the Taylor expansion 
Eq. H15|) becomes 



chem chem 

Per — n 0a 



chem 



2 ..chem 



dn„ 



-An n 



1 d 2 /4' 



2 dni 



-Ant 



(Al) 

In a degenerate non-magnetic system, the chemical po- 
tential (up to terms of order T 2 ) is given by 
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(67r 2 n CT ) 2 / 3 
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Truncating Eq. (|A1(I up to first order, is equivalent to 
request 
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which, using Eq. 1)A2|) becomes 
P{x) < 6. 



(A3) 



(A4) 



The above equation implies that, in the important case 
of a degenerate system, our approximation is valid up 
to injected polarizations of at least P(x) ~ 0.5^. In the 
special case of a single junction heterostructure it remains 
valid up to density polarizations of the order of 1 . This is 
shown in Fig.^J where the effects of the chemical poten- 
tial approximation discussed at the beginning of Sec. IV 
are presented for increasing electric field and FM polar- 
ization. Fig. 1131 compares the results obtained for R m in 
the FM/NMS structure using the non- interacting (degen- 
erate) approximation to the chemical potential, with the 
corresponding "exact" non-interacting degenerate results 
(see Appendix B). The figure shows that our approxi- 
mation for the chemical potential gives reasonably small 
errors (at most of the order of 10%) even for pf = 1. 
We caution though that this feature cannot be blindly 
generalized. 

Since |P(x)| < pf and the equality corresponds to the 
high field limit, we have analyzed systems in which pf 
satisfies the same limitations found for P(x). 
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FIG. 13: Correction to the (non-interacting) degenerate ap- 
proximation for R m , AR m /R m = (R% a - R% 9 ) / R% 9 vs elec- 
tric field for pf — 0.5 and Pf = 1 (as labelled). Here 
corresponds to the (non-interacting) "exact" degenerate cal- 
culation. 



eE-x + ei 



Fa 



(l±P(x)) 2/3 - 1 +5(B2) 



where the +(-) sign corresponds to a =| (j) and e Fa is 
the Fermi energy of the unperturbed system. 

Applying the conditions Q23 [l -H25 |) . we can now derive 
the equation for the injected density polarization P(0) 



- [1-P(0)] 2/3 } + ^ 1 pA T 1 ^=0. (B3) 



APPENDIX B: NON-INTERACTING 
DEGENERATE SYSTEM 

If we consider a non-interacting, degenerate NMS, 
Eq. (|18fl for the excess electrochemical potential becomes 

Afi a = eE ■ x + (e Fa - e° Fa ) + B (Bl) 



Once P(0) is known it is easy to derive all other quantities 
of interest. 

A similar scheme can be applied to the two junction 
system, where we obtain a system of two implicit equa- 
tions for the two quantities P(0) and P(xq), which must 
be solved self consistently . 
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